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ABSTRACT 

We report the detection of a z ~ 7 galaxy strongly lensed by the massive 
galaxy cluster Abell 2218 (z = 0.175) at 3.6 and 4.5 fim using the Spitzer Obser- 
vatory and at 1.1 fj,m using the Hubble Space Telescope. The new data indicate a 
refined photometric redshift in the range of 6.6 — 6.8 depending on the presence 
of Lya emission. The spectral energy distribution is consistent with having a 
significant Balmer break, suggesting that the galaxy is in the poststarburst stage 
with an age of at least ~ 50 Myr and quite possibly a few hundred Myr. This 
suggests the possibility that a mature stellar population is already in place at 
such a high redshift. Compared with typical Lyman break galaxies at z ~ 3 — 4, 
the stellar mass is an order of magnitude smaller (~ 1O 9 M ), but the specific 
star formation rate (star formation rate/M star ) is similarly large (> 10~ 9 yr _1 ), 
indicating equally vigorous star-forming activity. 

Subject headings: cosmology: observations — galaxies: formation — galaxies: 
evolution — galaxies: high-redshift — gravitational lensing 
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1. Introduction 

Locating and characterizing the first sub-galactic sources which may have been respon- 
sible for completing cosmic reionization and ending the "Dark Ages" represents the latest 
frontier in observational cosmology. Although gravitational instability theory allows the early 
formation of massive halos (e.g., Barkana & Loeb 2000), the complexities of gas cooling, star 
formation, feedback and the dumpiness of the intergalactic medium make observational pre- 
dictions highly uncertain (Stiavelli, Fall, & Panagia 2004). Data on representative early 
sources are needed to better understand when reionization occurred as well as what the 
luminosity and mass functions are of the sources responsible for the reionization. 

Searching for representative sources at redshift (z) of 6-10, the range considered to be 
the final stages of the reionization era, is a major observational challenge, not only because of 
the large luminosity distances involved but also since (as predicted by the hierarchical model 
of galaxy formation) the most likely systems to be found during this period have lower 
stellar masses than those at more moderate redshifts. The high magnification associated 
with strong gravitational lensing by massive foreground clusters offers a unique opportunity 
for detecting such galaxies. In recent years, this approach has been employed specifically to 
find galaxies during the reionization era, and has produced a number of successes (Ellis et 
al. 2001; Hu et al. 2002; Santos et al. 2004; Kneib et al. 2004a). The detection of a z ~ 10 
galaxy by Pello et al. (2004), however, is controversial (Weatherley et al. 2004; Bremer et al. 
2004). 

In this Letter, we return to analyze in more detail the physical characteristics of the 
triply lensed z ~ 7 source in Abell 2218 (z = 0.175) reported by Kneib et al. (2004a). The 
original redshift estimate was based on three pieces of observational evidence obtained for the 
two brighter components a and b: (1) reflection symmetry with respect to the well-defined 
z > 6 critical lines, (2) photometric redshift suggesting 6.6 < z < 7.1, and (3) a tentative 
spectroscopic detection of a continuum break at 9800 A, which, if due to Gunn-Peterson 
absorption at Lya, would indicate z = 7.05. Here, we report the Spitzer detection at 3.6 and 
4.5 /im of these two components as well as the new 1.1 fim detection by the Hubble Space 
Telescope (HST). With these new data, we will refine the photometric redshift and examine 
the physical properties of this lensed galaxy. 

2. Observations 

The 3.6 and 4.5 ^m images were obtained on UT 2004 January 2 with the InfraRed Array 
Camera (IRAC; Fazio et al. (2004)) on the Spitzer Observatory. In each filter, twelve 200- 



- 3- 



second images were taken with the small-step cycling dither pattern. The Basic Calibrated 
Data (BCD) images were combined using a custom IDL mosaicking routine with the final 
pixel scale of &!6 pixel -1 , half of the instrument pixel size. The IRAC observations also 
provided 5.8 and 8.0 /xm images of the same field simultaneously, but the z ~ 7 galaxy was 
not detected at these longer wavelengths. 

The 1.1 fj,m images were obtained on UT 2004 March 23 with the Near Infrared Camera 
and Multi-Object Spectrometer (NICMOS) on HST. A total of nine 2304 second (using the 
SPARS256 sequence) and one 1472 second (using the SPARS64 sequence) exposures were 
taken with the F110W filter using the NIC3 camera, which has a pixel scale of 0'.'2 pixel -1 . 

The central wavelengths and full widths at 20% of the transmission peak are listed in 
Table 1 for all the observed bands based on the total transmission curves. 



3. Results 

Figure 1 shows the ifST/NICMOS and Spitzer/IRAC images. Components a and b of 
the galaxy reported by Kneib et al. (2004a) are clearly detected at 1.1 /im (Figure la). At 3.6 
and 4.5 ^m, component b is clearly detected while component o is seen as a faint extension 
to the northwest of the z = 2.5 lensed submillimeter source SMM-A (SMM J16359+6612.6) 
discussed by Kneib et al. (2004b) (Figure lb and c). For comparison, the HST/NICMOS 
1.6 fim image by Kneib et al. (2004a) is also displayed (Figure Id), which shows the two 
components and SMM-A with a higher signal-to-noise ratio. When the light from SMM-A is 
subtracted using a two-dimensional elliptical Gaussian, component a shows up in the IRAC 
images (Figures le and f). 

At 1.1 /mi, photometry was performed in a manner consistent with that in Kneib et al. 
(2004a), giving the total fluxes of components a and b. In the IRAC bands, photometry was 
performed only for component b because the flux measurements with component a are highly 
uncertain due to the subtraction of SMM-A. For the photometry of component 6, we used 
an elliptical aperture with semi-major and semi-minor axes of 2" and 1", respectively, and 
aligned the major axis along the position angle of component b (23° East of North). The sky 
level was measured using similar ellipses to define the inner and outer boundaries of a sky 
area with a semi-major axis of 2 and 3.5 pixels, respectively. There is a slight sky background 
gradient along the east- west direction due to the bright source to the west, but the component 
b is far enough from this source that the gradient is quite linear across component b and 
the surrounding sky area. Therefore, we expect this gradient to be removed if we derive 
the local sky level as the average of all the sky pixel values. We also visually inspected the 
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image cross sections, and confirmed that the derived sky level is a good estimate around the 
position of component b. The photometric uncertainty is based on the scatter of sky pixel 
values, and our estimates are conservative in that the scatter not only includes the random 
sky pixel noise but is also inflated by the sky gradient along the east-west direction. 

Because of the large point spread functions at 3.6 and 4.5 /iin (FWHM ~ 1"7), the 
elliptical aperture misses a significant fraction of the total source flux. By convolving the 
F160W image with the IRAC PSFs, we determined that the flux measured with the elliptical 
aperture needs to be multiplied by a factor of 1.4 to account for this missing flux (i.e., aperture 
correction). 

The measured flux densities at 1.1, 3.6, and 4.5 ^m are listed in Table 1 together with 
the previous HST measurements presented in Kneib et al. (2004a). The 1.1 and 1.6yum 
images show a faint source north-west of component b. This source is also seen in the 
F606W, F814W, and F850LP images (Kneib et al. 2004a). The Spitzer images provide 
enough spatial resolution to conclude that there is very little flux emitted by this bluer and 
less distant source at 3.6 and 4.5 /im. 

4. Discussion 

The new photometric data permit us to construct the spectral energy distribution (SED) 
of component b, which we compare with the model SEDs simulated using the stellar popula- 
tion synthesis code GALAXEV (Bruzual & Chariot 2003). A grid of models was calculated 
with the parameters described in Table 2. The intrinsic (i.e., unmagnified) flux densities 
were assumed to be 25 times smaller than the observed values (Kneib et al. 2004a). For 
each model, the best-fit SED was sought by varying age and redshift and minimizing x 2 ■ 
The predicted fluxes were calculated from the model SEDs with appropriate total system 
throughput curves. The overall normalization of a model SED was determined by the F160W 
flux density measurement, which has the highest signal-to- noise ratio. This leaves five fil- 
ters (F814W, F850LP, F110W, IRAC 3.6 and 4.6^m) to constrain the SED shape. Given 
that little is known about the mode of star formation at z > 6, we tried to bracket various 
parameter ranges by calculating models for various star formation histories (see Table 2). 
The GALAXEV results were cross-checked against Starburst99 (Leitherer et al. 1999), which 
gave similar results. 

Figure 2a shows the best (i.e., minimum x 2 ) SEDs for each star formation history. The 
corresponding model parameters are listed in Table 2. All the model fits consistently give 
a redshift of 6.60-6.65. If we model the SED including Lya emission line with an observed 
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equivalent width of 100-300 A, the photometric redshift could increase up to z ~ 6.8. (This 
addition of a Lya line changes redshift, but does not affect the physical properties of the 
stellar population significantly.) The existence of such a Lya line is still consistent with the 
spectra presented by Kneib et al. (2004a) since there are a number of strong atmospheric OH 
lines in their wavelength range. A redshift above 6.8, however, is unlikely because it would 
require strong Lya emission, which would have been seen in the near-infrared spectrum in 
Kneib et al. (2004a). This suggests that the continuum break at 9800 A seen by Kneib et 
al. (2004a) is either spurious or an absorption line feature. 

Figure 2 also illustrates the importance of having Spitzer/IRAC photometry for such 
a high redshift galaxy. IRAC measures the SED above 4000 A in the restframe, thereby 
constraining the strength of the Balmer break. The IRAC measurements are indicative of a 
significant Balmer break (Figure 2), suggesting that the galaxy is observed well after the most 
vigorous stage of star formation (i.e., poststarburst), which is consistent with the fact that all 
the exponentially decaying star formation rate (SFR) models consistently give a galaxy age 
larger than the SFR e-folding time (r). A precise determination of the galaxy age, however, 
is difficult because it essentially depends on the star formation history we assume (Table 2). 
A reasonable lower limit on the age seems to be ~ 50 — 60 Myr based on the r = 10 Myr 
model (ignoring the instantaneous model, which is informative but unrealistic). The upper 
limit, on the other hand, is harder to constrain, and the 1 Gyr burst model is consistent 
with an age range of 300-700 Myrs (68% confidence interval), the large uncertainty being 
due to the slow SED evolution of this model in this age range. Taken as a whole, the model 
calculations in Table 2 suggest the existence of a mature stellar population whose age is at 
least 50 Myr and quite possibly a few hundred Myr. 

An interesting characteristic inferred from the model fits is an extremely blue restframe 
UV spectrum. The spectral slope index (3 (defined as fx oc X 13 ) is < —2, as is the case with the 
lower redshift UDF z-dropout galaxies for which NICMOS photometry is available (Stanway 
et al. 2004). Such a UV spectrum is consistent with a recent starburst and a normal IMF, 
but requires low extinction and/or low metallicity. Unfortunately, extinction and metallicity 
are degenerate, and are difficult to constrain independently with our photometric resolution. 
As an example, we show this degeneracy for the case of the r = 100 Myr model in Figure 2b 
with the corresponding model parameters in Table 2. 

By including the full range of models that produce acceptable fits in addition to those 
listed in Table 2, the various galaxy parameters are constrained to the following ranges 
(instantaneous models excluded): redshift: 6.6 — 6.8; age: 50-450 Myr; stellar mass: 0.5 — 
1 x 1O 9 M ; and SFR: 0.1-5 M yr _1 . These numbers reflect the ranges of the best-fit 
parameters without including the parameter uncertainties associated with each fit. The 
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mass is better constrained than the age, and the estimated mass (~ 1O 9 M ) is an order of 
magnitude smaller than those of typical Lyman break galaxies at z = 3 — 4 (~ 1O 1O M ) 
(Papovich, Dickinson, & Ferguson 2001; Shapley et al. 2001; Barmby et al. 2004). However, 
the specific star formation rate (SFR/M star ) (e.g., Brinchmann et al. 2004) is similarly large 
(> 10 9 yr _1 ), indicating equally vigorous star- forming activity. 

The possible detection of a mature stellar population in this galaxy allows us to extend 
the argument made by Kneib et al. (2004a) on the basis of the high surface density on the 
sky of sources similar to that studied here. As implied by its discovery in a small survey 
area, unless we have been extraordinarily fortunate, such sources are likely to have a mean 
surface density of ~ 1 arcmin -2 and remain luminous over an extended period of cosmic 
history. Accordingly, it is likely that they contribute significantly to the UV photon budget 
necessary for cosmic reionization. 

5. Conclusions 

The new Spitzer and HST data indicate that the redshift of this lensed galaxy is indeed 
very high, probably at z ~ 6.6 — 6.8. Comparison with a variety of stellar population 
synthesis models indicates that the galaxy is in the poststarburst stage with an age of at 
least ~50 Myr and quite possibly a few hundred Myr, which suggests the possibility that a 
mature stellar population is already in place at such a high redshift. Unless we have been 
fortunate in its discovery, it is likely that such sources with extended lifetimes contribute 
significantly to cosmic reionization. 

We acknowledge helpful discussions with G. Smith, D. Stern, R. Pello, D. Schaerer, and 
thank G. Neugebauer for commenting on the manuscript. JPK acknowledges support from 
Caltech and CNRS. The study of Abell 2218 as a cosmic lens is supported by NASA STScI 
grant HST-GO-09452.01-A. 
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Fig. 1. — HST/N1CM0S and Spitzer/IRAC images of the z ~ 7 lensed pair: (a) New 
NICMOS 1.1 fim image. Components a and b as well as the z=2.5 submillimeter source 
SMM-A are marked. The dashed line indicates the z > 6.5 critical curves; (b) IRAC 3.6 
jum image; (c) IRAC 4.5 /im image; (d) NICMOS 1.6 jum image presented in Kneib et al. 
(2004a); (e) IRAC 3.6 ^m image with SMM-A subtracted; (f) IRAC 4.5 /im image with 
SMM-A subtracted. Each image is 16" on a side. In the panels (a) and (d), north is up 
and east is to the left. In the panels (b), (c), (e), and (f), we preserved the original array 
orientation, which is 11° rotated counter-clockwise from that of the panels (a) and (b), in 
order to display the results of the SMM-A subtraction accurately. The dashed-line ellipses 
shown in (b) and (c) around component b indicate the areas used for photometry. Similar 
ellipses are drawn around component a to mark its position. 
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Fig. 2. — SED model fits to the observed SED of the z ~ 7 galaxy (component 6): (a) The 
best (i.e., minimum model for each star formation history; (b) the exponentially decaying 
starburst model (r = 100 Myr) for a range of metallicities. The restframe wavelength at 
z = 6.65 is also shown. The error bars indicate the AA and flux density uncertainty of each 
band shown in Table 1. The measured flux at 1.1 jttm is expected to be significantly lower 
than the true continuum level because the F110W filter passband extends below 1216 A in 
the restframe. 
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Table 1. Flux densities of components a and b 



Filter A c a AA (20%) a Flux density 

a b 





(H 


(H 




(^Jy) 


F8UW 


0.801 


0.242 


0.09±0.02 


0.03±0.01 


F850LP 


0.915 


0.214 


0.39±0.02 


0.30±0.02 


F110W 


1.121 


0.604 


0.78±0.05 


0.67±0.04 


F160W 


1.612 


0.418 


0.87±0.04 


0.80±0.04 


3.6 jttm 


3.561 


0.773 




1.2±0.3 


4.5 fim 


4.510 


1.015 




1.0±0.2 



a The effective wavelength and the full width at 20% of 
the peak calculated from the total transmission curves. 



Note. — The IRAC flux densities of component a are 
highly uncertain due to the subtraction of SMM-A, and 
therefore are not listed. 



Table 2. Model parameters 



Model 




Metallicity 


z 


Age a 


Mass 


TV 


SFR 


Z form 






(Z©) 




(Myr) 


(10 s M ) 




(M /yr) 








Minimum 


X models 










Instantaneous burst 




0.4 


6.60 


35 (28-40) 


5.5 


0.25 




6.8 


e-decaying burst (r = 


10 Myr) 


0.02 


6.60 


64 (56-72) 


8.2 


0.50 


0.2 


7.0 


e-decaying burst (r = 


50 Myr) 


0.4 


6.60 


102 (82-114) 


6.3 


0.25 


2.2 


7.3 


e-decaying burst (r = 


100 Myr) 


1.0 


6.65 


181 (150-201) 


8.3 


0.0 


2.0 


8.1 


1 Gyr burst 




1.0 


6.65 


453 (287-693) 


9.7 


0.0 


2.6 


12.3 




e-dccayin 


g (r = 100 Myr) models with sub-solar metallicities 






e-decaying burst (r = 


100 Myr) 


0.005 


6.60 


161 (121-185) 


11 


0.75 


3.3 


7.8 






0.02 


6.60 


143 (124-175) 


8.5 


0.50 


3.2 


7.7 






0.2 


6.60 


90 (66-120) 


5.4 


0.50 


4.2 


7.2 






0.4 


6.60 


128 (105-156) 


6.2 


0.25 


2.8 


7.5 



a The numbers in parentheses indicate the 68% confidence interval based on A% 2 . 



Note. — Models with GALAXEV (Bruzual & Chariot 2003) explored the following range of parameters: (1) 
Initial mass function (IMF): Salpctcr with the lower and upper mass cutoffs of 0.1 and 100 Mq; (2) metallicity: 
0.005, 0.02, 0.2, 0.4, 1.0, and 2.5 Z Q ; (3) dust extinction: t v = 0.0, 0.25, 0.5, 0.75, 1.0, 1.25, and 1.5, where t v is 
the total effective V-band optical depth seen by young (age < 10 7 yrs) stars. The optical depth toward old stars 
(age > 10 7 yrs) was set to be 1/3 of this. See Chariot & Fall (2000) and Bruzual & Chariot (2003) for more detail; 
and (4) star formation histories: instantaneous burst, exponentially decaying burst (SFR oc e~'/ T , where t is the 
galaxy age and t is the SFR e-folding time scale) with r = 10, 50, and 100 Myr, and 1 Gyr burst (a constant SFR 
model with the SFR scaled such that all the gas gets consumed in 1 Gyr). The flux below Lya (A res t = 1216 A) 
was set to zero, and no Lya emission line was included. SFRs listed are intrinsic model values. We assumed a 
magnification factor of 25 (Kneib et al. 2004a) and the cosmological parameters of f!in = 0.3, Qa = 0.7 and Ho=70 
km s _1 Mpc -1 . 



